NASA Contractor Report 4201 


Gear Optimization 


G. N. Vanderplaats, Xiarig Chen, 
and Ning-tian Zhang 


GRANT NAG3-683 
DECEMBER 1988 


NASA 


NASA Contractor Report 4201 


Gear Optimization 


G. N. Vanderplaats, Xiang Chen, 
and Ning-tian Zhang 
University of California 
Santa Barbara, California 


Prepared for 

Lewis Research Center 

under Grant NAG3-683 


f\l/\SA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Division 


1988 



CONTENTS 


Page 


LIST OF SYMBOLS iv 

1.0 INTRODUCTION 1 

2.0 NUMERICAL OPTIMIZATION CONCEPTS 2 

3.0 THE COPES/ADS OPTIMIZATION PROGRAM 4 

4.0 SPUR GEAR ANALYSIS AND OPTIMIZATION 5 

4. 1 THE SPUR GEAR PROGRAM 20 

4.2 PARAMETERS CONTAINED IN THE GLOBAL COMMON BLOCK 2 1 

4.3 DESIGN EXAMPLES 27 

5.0 SPIRAL BEVEL GEAR ANALYSIS AND OPTIMIZATION 31 

5 . 1 THE SPIRAL BEVEL GEAR PROGRAM 44 

5.2 PARAMETERS CONTAINED IN THE GLOBAL COMMON BLOCK 44 

5.3 DESIGN EXAMPLES 48 

6.0 SUMMARY 53 

7.0 REFERENCES 55 


preceding Page 


BtANK NOT FILMED 


III 


mm 



Symbols 


E l’ E 2 

^2 


V V 
v r 2 


normal load (lb) 

pressure angle (deg) 

tangential load (lb) = W n • costj) 

the number of pinion teeth 

the number of gear teeth 

diametral pitch (1/in) 

base pitch, (in/tooth), = 7tcos<t>/P 

pitch radius of pinion (in) r^ = Nj/2P 

pitch radius of gear (in) r 2 = N^P 

addendum circle radius of pinion (in) r j = ^ + 1/P 

addendum circle radius of gear (in) = r 2 + 1/P 

base circle radius of pinion (in) r fel = ^ • cos§ 

base circle radius of gear (in) r fe2 = r 2 • co$4> 

center distance (in) C = ^ + r 2 

addendum (in) a = 1/P 
dedendum (in) b = 1.25/P 
face width of pinion (in) 

face width of gear (in) 

face width in contact (in) f = min(f 1 ,f 2 ) 

density of pinion (lbm/in ) 

density of gear (lbm/in ) 

Young's modulus 
Poisson’s ratio 

torsional modulus of elasticity 
speed of pinion (rpm) 

speed of gear (rpm) n g = Nj/N 2 • n p 
pitch line velocity (in/sec) V = 71^/30 • n p 
curvature radius (in) 

rolling velocity at point of contact (in/sec) 

Vj = itn p /30 ‘Pj V 2 = TCn g/30 • p 2 
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1.0 INTRODUCTION 


The design of gears is a highly complicated task, and the need to develop light 
weight, quiet and more reliable designs has resulted in a variety of changes in the design 
process. Very little published work is available on the subject of gear optimization. One 
exception is the work of Savage and co-workers [1], who present a method for the 
optimum design of gears with the object of minimizing size and weight. Their approach 
consists of comparing several candidate designs to obtain the best design and clearly this 
provides little assurance of obtaining a true optimum design. Also, identifying the best 
design using graphical methods becomes tedious once the number of design variables 
exceeds two. Therefore, a more general and systematic approach to gear design is 
desirable. 

Here, mathematical programming techniques, commonly referred to as numerical 
optimization methods, are investigated to provide a reliable design methodology for gears. 
These techniques offer a logical approach to design automation and also can handle a wide 
variety of design variables and constraints which are difficult to visualize using graphical 
methods [2]. 

As part of this study, two gear optimization programs have been developed, one for 
spur gears and one for spiral bevel gears. The programs evaluate a wide variety of 
functions required by the optimizer. Typical design variables include the number of teeth, 
diametral pitch and face width. Design objectives may be minimum weight, minimum 
dynamic load factor or maximum gear life, as examples. Constraints may be imposed on 
loading or Hertzian stress, volume, and tooth velocity as examples. Also, it should be 
noted that objective and constraint functions are interchangeable. It is possible to minimize 
dynamic load factor with a required life or, alternatively, to maximize life with a bound on 
load factor. 

In this report, the basic concepts of numerical optimization, as a gear design tool, are 
outlined first. Next, the COPES/ADS optimization program is briefly described. This 
program was used as the optimization tool throughout this study. The particular advantage 
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of this program is that the user is free to choose different design variables, objective 
function and constraint functions using input data when the program is executed, without 
modifying the program itself. The analysis capabilities developed here for spur gear and 
spiral bevel gear design are described next. It is important to note that these programs 
perform analysis only, and can be used for that purpose, without optimization. They have 
been written in the form needed to couple them to the COPES/ADS program when 
optimization is desired. When discussing each of these analyses capabilities, design 
examples are given to demonstrate their application in the optimization mode. Finally, the 
results of this project are summarized and continuing research efforts are described 

2.0 NUMERICAL OPTIMIZATION CONCEPTS 
The constrained optimization problem is mathematically stated as follows [1] 


Minimize (or maximize) 

F(X) 


(1) 

Subject to 

gj(X)£0 

j = l,m 

(2) 


h k (X) = ° 

k = 1,/ 

(3) 


xf^Xi^x" 

i = l,n 

(4) 


Where, F(X) is referred to as the objective function. gj(X) are inequality constraint 
functions and h k (X) are equality constraint functions. X., i=l,n, are design variables. 
Equation (4) defines side constraints which impose limits on the region of search for the 
optimum. 

A typical optimization problem for spur gear design can be stated as 


Maximize 
Subject to 
limits on 

I 

I 


Life 

1) Dynamic load 

2) Bending strength 
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3) Surface durability 

4) Scoring 

5) Weight 

6) Gear size 

7) Contact ratio 

8) Involute interference 

9) Width to diameter ratio 

Alternatively, the dynamic load factor might be minimized with a lower bound 
constraint on life. Using COPES [3] and ADS [4,5], only input data needs to be changed 
for different problems. 

Design Variables 

The key element in formulating an optimization problem is the selection of the 
independent design variables that are necessary to characterize the design of the system. 
Normally, it is good to choose those variables that have a significant impact on the 
objective function. Based on this criterion, the number of teeth on the pinion and gear, Nj 
and N 2> the diametral pitch P and facewidth of the pinion and gear, f } and f 2 are typically 
chosen as the design variables. Since present optimization techniques do not deal well with 
integer variables, and N 2 are treated as a continuous variables in the optimization. The 
final answer can be rounded to the nearest integer value. 

In principal, any input parameter to the gear analysis program can be taken as a 
design variable. The input data to the COPES/ADS program defines the number and 
identity of the design variables at run time. Also, design variables may be "linked" so that, 
for exampe N 2 = 3 * Nj would assure a 3 to 1 gear reduction. 

Objective Functions 

The surface pitting fatigue life of a gear set in mesh and the dynamic load on a gear 
tooth are the objective functions in life maximization problem and dynamic load 
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minimization problem, respectively. Other objective functions may be size, weight, or 
tooth velocity as examples. 

Constraints 

Just as any response quantity calculated by the gear analysis program can be treated 
as the design objective, so can any quantity by treated as a constraint. As listed above, 
typical constraints here include limits on bending strength, surface durability, scoring, 
weight, size, contact ratio, interference and width-to-diameter ratio. Also, life and dynamic 
load can be treated interchangeably as objective or constraint functions. All that is needed 
to change the definition of the objective or constraint functions or the bounds on the 
constraints, is to change the input data to the COPES/ADS program. 

3.0 THE COPES/ADS OPTIMIZATION PROGRAM 

The COPES/ADS optimization program [3-5] was used in this study as the design 
tool, coupled with the spur gear or bevel gear analysis. The ADS program is a general 
purpose optimization program containing a variety of modern nonlinear programming 
algorithms. The COPES program is a control program which acts as the interface between 
ADS and the user-supplied analysis. 

The basic concept of using COPES/ADS is that the user must supply an analysis 
program in subroutine form, called ANALIZ. The analysis is separated into three parts, 
being INPUT, EXECUTION, and OUTPUT. All parameters in the analysis that may be 
design variables, objective functions or constraints are contained in a single labeled 
common block called GLOBCM. The ordering of the information in GLOBCM is 
arbitrary. 

The entries into the GLOBCM common block are treated as a set of sequential 
information to be accessed by COPES and the input data to COPES refers to these entries 
by their location. For example, the input to COPES may tell the program to minimize the 
value of entry number 3, while locations 1, 6 and 21 define three design variables. Finally, 
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locations 30 through 35 and locations 7 and 2 can be constrained. The details of how the 
data is input to the program are given in reference 3. 

The key to using the COPES/ADS program is that it provides the designer a great deal 
of flexibility in choosing the design variables, objective and constraints, and that no 
program modifications are needed between design runs. It is only necessary to change the 
input data. This allows the designer to compare a variety of optimized designs, as opposed 
to the usual single variable trade-off studies. Also, the program is in no way limited to the 
analysis of spur and bevel gears described here, but is applicable to any design problem for 
which to ANALIZ subroutine is avialable. 

4.0 SPUR GEAR ANALYSIS AND OPTIMIZATION 

One of the design capabilities developed as part of this study is for spur gears. The 
study began by converting the TELSGE program, provided by NASA Lewis Research 
Center, to a subroutine ANALIZ and gaining experience in design for life maximization. 
This gear program does not include many of the constraints that must be considered by 
optimization, such as bending strength and scoring. However, TELSGE provided an 
excellent starting point for development of a more detailed program for use in optimization. 
The basic analysis capabilities of that program are described here and examples are given to 
demonstrate its capabilities. The purpose here was not to develop new analytical models 
for gears, but rather to identify and program good models into a general design program. 
The program is written in modular form so that any of the component parts can be modified 
and updated in the future. 

Analytical Model for the Life of Spur Gears 

The life of spur gears is calculated from an analytical model developed by NASA 
Lewis Research Center, which is programmed in subroutine LIFESG [6-8]. The basic 
gear geometry with geometrical parameters is shown in figures 1 through 4 . 
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Here, it is assumed that the life of spur gears is mainly determined by the surface 
pitting fatigue failure caused by repeated applications of high surface contact stress. 

The life of a single pinion tooth in millions of pinion rotations with 90 percent 
probability of survival is given as 





Csin<j> 

r bl+Ll (Csio * -r bl'*Ll ) 


( 6 ) 


6 L1 is the roll angle corresponding to the lowest point of single tooth contact on the pinion 
where the maximum Hertzian stress is considered to occur 



In equation (5), l { is the involute profile arc length of the pinion (in) 

l l = r bl ^H1 ( 5 1 + Pli + 2 0hi^ (10) 
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where, 


0H.= 


2p, 


- <V r al - r bl + V r a2 - r b2 - C si" 0» 


bl 


(ID 


Here, ^ is calculated as the involute length across the heavy load zone as a justifiable 
simplication. 

For a single gear tooth, the life in millions of pinion rotations with 90 percent 
probability of survival 


L 2p= ( ¥7 )( i r ' L > 


( 12 ) 


where, e is the Weibull exponent which is taken as 2.5. I 2 is the involute profile arc length 
of the gear (in) which is calculated similarly to ^ 


l 2 r b2 P H2 ^2 + + 2^ H2^ 


(13) 


5„ = 


C sin * - J - r 


b2 


(14) 


_ 7 r al ~ r bl + V f a2 ~ r b2 ~ C sin ♦ ~ P b 

Pl2 r b2 

R 2 Pb ~ ( 7 r al ~ r bl + V r a2 - r b2 ~ C sin ♦> 
Ph2 " r b2 


(15) 


(16) 
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For a pinion alone, the life in millions of pinion rotations with 90 percent probability of 
survival 


L 


P 



(17) 


For a gear alone. 


L 


Gp 




(18) 


The life of gear set in mesh in 
survival is 

or 


millions of pinion rotations, with 90 pecent probability of 




L M = f N i<r-) e+N 2<L7 

1 2p 


) e ) 


(19) 


Dynamic Load 

Dynamic load on a spur gear tooth is calculated by using Buckingham’s solution 
[9,10], expressed as 


W <r W . + V f a < 2f 2 - f a> 


( 20 ) 


where, f is the acceleration load which is a function of f . , the force needed to accelerate 
rigid bodies through the critical error, and f 2> the force required to do the work of 
deforming the tooth by elastic deformation equal to the amount of the critical mesh error, 

f a = f l* f / f . +f 2> (21) 
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( 22 ) 


fj = [tan<t> (1 - cos<|))/l 50<J> 2 ] [(l/r^ + (l/r 2 )] mV 2 = HmV 2 

f 2 = W t [(e/d) + l] (23) 

and m is the effective mass at pitch line of gears 
e is the measured error on pair of mating teeth 
d is the total deformation of mating teeth under applied load 

For 20-deg full-depth form, 


d = 9. 00 (W /f ) [-T7- + -=r-] 
a l *2 


(24) 


In equation (22) 


m = nijiryOnj + m 2 ) (25) 

m l = m p + n^ (26) 

m_ = m + m. (27) 

2 g d 

Where, nij and m 2 are effective masses acting at pitch line of pinion and gear, respectively. 

nip and m g are effective masses of pinion blank at r 1 and gear blank at r 2 , 
respectively. 


m 


P 


2 ' 



(28) 


9 


K B 4 , 

m g = — r a 2 - f 2 

m b = ( a /b^ + 4A 1 C 1 - Bj) /2Aj 
A = Hm V 2 

1 a i 

B J = (m p + m 2 ) A J + em 2 Z 1 
c l = « n, . 1 - m 2 Z l 

m d = <7 B 2 +4A 2 C 2 - B 2> /2A 2 

A- - Hm. V 2 
2 a 2 

B 2 = (rn g + A 2 + em^ 

C 2 = em a 2 m i Z 2 
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( 29 ) 


(30) 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

ar Tj and r 2 , respectively. 

(38) 


( 39 ) 


where, d f and d f are diameters of shaft corresponding to pinion and gear, respectively. 

r i *2 

Lj and L 2 are length of shaft corresponding to pinion and gear, respectively. 

Zj and Z 2 are elasticity factors of shaft corresponding to pinion and gear, 
respectively. 


Z l = 7td 4 f • G/(32r*- LJ 


(40) 


Z 2 = * d f 


G/(32r2- L 2 


(41) 


Subroutine DYLOAD is programmed to calculate the dynamic load. 


Bending Strength 

In order to avoid gear tooth breakage, the bending stress of a gear tooth may not 
exceed the allowable tooth stress. 

The bending stress is calculated by using the AGMA formula [11]. 


W P 

CT= K^T 


(42) 


where, K y is dynamic factor. 

J is AGMA geometry factor. 

Since the dynamic load W d can be calculated by equation (20), equation (42) can be 
converted to 


W P 
a 

~ fj 


(43) 
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AGMA geometry factor J for external spur gear teeth is determined from the analytical 
model presented by R.G. Mitchiner and H.H. Mabie [12]. 


J = 


Y 

K f - m N 


(44) 


where, 

Y is the Lewis form factor. 

K f is the stress concentration factor, 
rtij^ is the load-sharing ratio. 


The rr^ may be taken as the profile contact ratio mp 


m N = m ] 



K COS <|> 


- c sin <|>)P 


(45) 


Kf is defined as 


where, 


Kf =H + (T) L + ( T 

c y e 

(46) 

H = 0.34 - 0.4583662 $ 

(47) 

L = 0.316 - 9.4583662 <> 

(48) 

M = 0.290 + 0.4583662 <|> 

(49) 

r = r f +(b-r f ) 2 /(R+b-r f ) 

(50) 


R = pitch radius. 
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r f = fillet radius on hob tip. 

T = tooth width between points of tangency of parabola and root fillets. 

R = radius from center of gear to apex of parabola, 
c 


1 7t 

r , = : — r (~j — cos d> - b sin A) 

f 1 - sin <j> 4p Y T 


( 51 ) 


R c = 


R cos <j> 


cos(<() A - [~ + inv(J) - inv<t> A ] ) 


(52) 


inv<b = tan<|> - f 


(53) 


cosQ] 2 

. . _i R+ a N 

^ A = tan ( R ~ > 

rTT cos0 


(54) 


where, 0 is calculated by using Newton's method. 


0 = = e. - 


fOj) 


i+l i f(0 .) 


(55) 


f(0) m j£ + X 


2( R c -y E > 


(56) 


b-r f 

^(-^tantf + O) 

m. = r 

JE b-r f 

("^q ) - tan (P + 0) 


(57) 
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(58) 


0 = 


TC_ 

N 


X E = (R - b + r f ) sin (3 + 0) - R0 cos (3 + 0) - 
[(b - r f )sin(P + 0) + R0 cos(P + 0)] 
y E = (R-b+ r f ) cos (3 + 0) + R0 sin (3 + 0) + 
[R0 sin(P + 0) - (b - r f ) cos (3 + 0)] 


J~( b-r f ) 2 + R 2 0 2 


J(b- T[ ) 2 + R 2 e 2 


(59) 


(60) 


f (0) is calculated by using first-forward differential method. In equation (46), 


T = 2X 

E 


(61) 


Y is defined as 


Y = 


P 

°”*L .l5 _ “V 
cos (}> 1 X T 


(62) 


where, 


X = 



R 


c -y 


E 


(63) 
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( 64 ) 


♦ 


/ i-r-R 
4 -lW l R + 

L = tan ( g- 


COS <()] 


R + a 


a 

cos <}) 


- inv[tan 


fl ITS: 

-1,V l R + ; 


COS (J)]' 


-)] 


R + a 


cos ()) 


+ inv<)> . 

AGMA geometry factor J is calculated in subroutine GEOFAC. The bending stresses 
of pinion and gear tooth are calculated in subroutine BENDST. 

The constraint condition for bending stress is 


a £ a . 


where, a is material allowable stress. 


Surface Durability 

The gears must be designed so that the dynamic surface stresses are within the 
surface endurance limit of the material. 

The surface contact stress (Hertzian stress) [11] 


a 


H 



w d 

Ttf COS <)) 


q/ Pl ) + (i/p 2 ) 

[<1 - M-f ) / Ej] + [(1 - M- 2 ) / E 2 ] 


(65) 


where, and p 2 are calculated at pitch point. The surface fatigue strength for steels is 


S c = (0.4 H b - 10) x 10 3 psi 


( 66 ) 


H d is the Brinell hardness of the softer of the two contacting surfaces. 
£> 

Modified by effective factors, the corrected fatigue strength 
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( 67 ) 


C C 

c _ L^H q 

H “ C T C n C 
T R 


where, = life factor 

C H = hardness - ratio factor 
Cj. = temperature factor 
= reliability factor. 

The constraint condition for surface durability is 


°H“ S H • 


C„ and S„ are calculated in subroutine SURDUR. 

H ri 

Scoring 

There are three constraint conditions to be considered for scoring failure, 
flash temperature limit, scoring criterion number and specific film thickness. 

The flash temperature is expressed as [13,14] 


T = T 

if i 


+ __£fW 

» cos * • 


where, T f = flash temperature, *F. 

T fe = temperature of blank surface in contact zone, *F. 

Cf = material constant for conductivity, density and specific heat. C f 
f f = coefficient of friction, f^ = 0.06. 

B = width of band of contact. 


They are 


( 68 ) 


0.0528. 
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(69) 


B = 0. 00054 [ 


W tPlP 2 

cos <|> • f (pj +p 2 ) J 


0. 5 


where, p ( and p 2 are calculated at the contact point where the flash temperature is higher. 


The scoring criterion number = 


W t 3/4 

(r-f-y 



(70) 


The critical scoring-criterion number 


= 60T fe + 23000. Which is formulated from the table [14]. (71) 

The specific film thickness [13] 


X = 


min 

a 


(72) 


o = Jo* + V\ 


(73) 


where, h . = minimum oil film thickness, in. 
mm 

Op = pinion average roughness, rms. 

O = gear average roughness, rms. 

8 

The minimum elastohydrodynamic film thickness is [13,15] 


h . = 
mm 


2 65a a54 (v n u) a7 E’ a03 R ft43 


W 


0.13 


(74) 
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where, 


a = pressure viscosity coefficient. 
v q = absolute viscosity. 


E' 


1 l -» 2 l l -*l 

2 ^ E + c * 


( 75 ) 


1 


U= 2 (U 1 + U 2 ) 


(76) 


Uj and u 2 are rolling velocities of pinion and gear at point of contact, respectively. 


Ttn. 


u l " 30 ‘ p l 


(77) 


7m , 


u 2 ~ 30 * p 2 


(78) 


R = 


P,P 


1 K 2 


Pl + P 2 


(Pj and p 2 are calculated at the lowest point where the worst Hertzian stress is 
considered to occur.) 


W t 

W = specific loading. It is taken as * here. 

To avoid the scoring failure, the flash temperature should be below a certain limit, the 
specific film thickness may not be less than 1.0, and the scoring number may not exceed 
the critical number. 

These functions are calculated in subroutine FLASH, SCORING, FILMTH, 


18 


FLMFC. 


Gear Weight and Size 

In order to increase the life of gears or reduce the dynamic load on gear tooth without 
the cost of increasing the weight and the gear size, constraints may be placed on the weight 
of the pinion, gear and the gear mesh center distance C. Alternatively, one of these 
parameters may be taken as the design objective. 


2 

Weight of the pinion = p A • Jtr l • f 1 

(79) 

Weight of the gear = P B r 2* f 2 

(80) 

Gear mesh center distance, c = ^ + r 2 

(81) 


Contact Ratio 

The contact ratio 

m P n cos <{> (82) 


For standard spur gears, m 


should be greater than 1.4 and is normally less than 2.0 


[ 11 , 12 ]. 


Involute Interference 

Involute interference occurs when the driving gear tooth makes contact with the 
driven tooth before the involute portion of the driving tooth comes within range. That is, 
when contact occurs below the base circle of the driver gears on the non-involute portion of 
the flank. 
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For no involute interference, a constraint condition is imposed [1]. 



Width to Diameter Ratio 

For the tooth load to be uniform, a common criterion is to limit the length to diameter 
ratio of the line of contact [1]. 

This ratio is 


X R = ^ (85) 

An additional requirement is Nj <; N 2 , since Nj and N 2 are corresponding to pinion 
and gear, respectively. 

4 . 1 The Spur Gear Program 

The detail structure of COPES and ADS will not be discussed here, since it has been 
already presented in References 3-5. The spur gear analysis subroutines include ANALIZ, 
TREAT, BENDST, GEOFAC, GEOFUN, GEODF, DYLOAD, SURDUR, SCORIN, 
FLASH, FILMS C, WEIGH, RATIOC, LIFESG, FILMTH and OUTPUT. ANALIZ is 
the subroutine required by COPES such that the analysis subroutines could be coupled to 
COPES. TREAT evaluates required items in terms of the initial values of design variables 
and constant data input by users. BENDST evaluates the bending stresses by using 
AGMA formula, where the required geometry factor J is calculated in routine GEOFAC. 
GEOFUN and GEODF evaluate functions and function gradients required in geometry 
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factor calculations, respectively. DYLOAD evaluates the dynamic loading on gear teeth 
using Buckingham's solution. SURDUR calculates the Hertzian stress to observe surface 
durability, while the scoring criterion number and the flash temperature are determined in 
SCORIN and FLASH, respectively. The ratio of film thickness to composite surface 
roughness is calculated in FILMSC, which is one of the parameters in scoring constraint 
conditions. WEIGH calculates the weight of gear mesh and RATIOC calculates the contact 
ratio. LIFESG calculates surface fatigue life of spur gear mesh, which is one of the 
objective functions. FILMTH calculates the minimum elastohydrodynamic film thickness. 
Finally, OUTPUT prints the results at each step. 

4.2 Parameters Contained in the Global Common Block 

In order to perform design operations, the COPES program must access the data in 
common block GLOBCM. The location in GLOBCM where a specified parameter resides 
must be defined in the COPES input data. The statement of the common block GLOBCM 
in this spur gear optimization design program is 

COMMON /GLOBCM/ LFM,WD,NTA,NTB,PD,FWAD,FWBD,CTAA, 

DSTBD(2),SST,CSORN,TEMFL,LAMDA, RATIO, WEIGHT(2),CENL, 

RATLD(2),AINTF,DNT 

Where, LFM is the life of gear mesh in millions of pinion rotations, which is the value 
of the objective function in the life maximization problem. 

WD is the dynamic load on gear tooth (lb.), which is the value of the objective 
function in dynamic load minimization problem. 

NTA and NTB are the number of teeth of pinion and gear (in.), respectively, which 
are design variables. 

PD is the diametral pitch (1/in.), which is one of the design variables. 

FWAD and FWBD are the face width of pinion and gear (in.), respectively, which 
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are design variables. 

CTAA is the pressure angle (deg.), which is an input constant. 

DSTBD is an array having two elements. DSTBD(l) and DSTBD(2) are the 
maximum bending stresses acting on the tooth of pinion and gear (psi.), respectively. 
Here, the dynamic load is considered in the calculation. 

SST is defined as the surface contact stress (Hertzian stress) minus the surface fatigue 
strength (psi.). 

CSCORN is defined as the scoring criterion number minus the critical scoring 
criterion number. 

TEMFL is the flash temperature (*F). 

LAMDA is the specific film thickness. 

RATIO is the contact ratio. 

WEIGHT is an array having two elements. WEIGHT(l) and WEIGHT(2) are the 
weight of pinion and gear (lb.), respectively. 

CENL is the gear mesh center distance (in.), which presents the size of the gear 
mesh. 

RATLD is an array having two elements. RATLD(l) and RATLD(2) are the ratio of 
width to diameter of pinion and gear respectively. 

AINTF is a constraint condition corresponding to involute interference, which is 
defined as 

AINTF = c sin 6 - [r\ — r£ 

V *2 b 2 


AINTF should be greater than 0, since 


c sin <{> 


> 




2 
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DNT is defined as NTB minus NTA. DNT should not be less than 0, since the gear 
is always not less than the pinion. 

DSTB, SST, CSCORN, TEMFL, LAMDA, RATIO, WEIGHT, CENL, RATLD, 
AINTF and DNT are values of constraint functions. The definitions of some terms have 
been explained previously. 

Input Data Description 

There are two parts of input data file. The first part corresponds to COPES; and the 
second part corresponds to the spur gear analysis subroutines. The details of the COPES 
input may be found in the COPES manual (3) and so are not described here. 

Spur Gear Analysis Input Data 

There are six rows of data in the second part input file. The form of the input data file 
is as follows; 

NTA,NTB,PD,FWAD,FWBD,CTAA 

TTLD,CYCLP 

DENSA,DENSB,YMAD,YMBD,EPSA,EPSB,SRMS1,SRMS2 
ERR.SHL1 ,SHL2,SHD 1 ,SHD2 
HARDB,CL,CH,CT,CK 
ALPHA, VISD, BETA, TEMLUB 

Where, NTA and NTB are the number of pinion and gear teeth, respectively. 

PD is the diametral pitch (in. ~ 

FWAD and FWBD are the face widths of pinion and gear, respectively (in.). 

CTAA is the pressure angle (deg.). 

TTLD is the normal load acting on gear tooth (lb.). 

CYCLP is the pinion speed (r/s). 

DENSA and DENSB are the material density of pinion and gear, respectively 
(lbm/in 3 ). 
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YMAD and YMBD are Young's modulus of pinion and gear, respectively (psi.). 
EPSA and EPSB are Poisson's ratio of pinion and gear, respectively. 

SRMS1 and SRMS2 are surface roughness corresponding to pinion and hear, 
respectively (microin.) 

ERR is the measured error on pair of mating teeth (in.). 

SHL1 and SHL2 are the lengths of shaft corresponding to pinion and gear, 
respectively (in.). 

SHD1 and SHD2 are the diameters of shaft corresponding to pinion and gear, 
respectively (in.). 

HARDB is the Brinell hardness, which is required in scoring observation. 

CL, CH, CT, CK are the life factor, the hardness ratio factor, the temperature factor 
and the reliability factor, respectively. These factors are required in surface durability 
observation. 

ALPHA is the pressure viscosity coefficient. 

VISD is the viscosity (lb. s/in ) at ambient temperature. 

BETA is the temperature viscosity coefficient. 

TEMLUB is the ambient temperature (deg. F). 

ALPHA, VISD, BETA and TEMLUB are required input data in film thickness 
calculation. 

All of above input data are real number, including the number of teeth NTA and 

NTB. 

Input Data Example: Life Maximization 

The input data for a design example of the life maximization problem is presented in 
the following Table. 
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Table 1 


$ BLOCK A: 

OPTIMIZATION DESIGN FOR LIFE OF SPUR GEARS 
$ BLOCK B: 

2.3 

$ BLOCK C: 

0,5,7,1000 
$ BLOCK D: 

$ BLOCK E: 

$ BLOCK F: 

5,1, +1.0 
$ BLOCK G: 

30.. 50.. 36. 

8 .. 12 .. 8 . 

0.2,0.4,0.25 

$ BLOCK H: 

1.3 

1.4 

2.5 

3.6 

3.7 

$ BLOCK I: 

10 

$ BLOCK J: 

9.10.1 

0.,0., 180000., 0. 

11 . 11.1 

0.,0.,1.E15,0. 

12,12,1 

-l.E+10,0.,0., 

13.13.1 
0.,0.,275.,0. 

14.14.1 

1.. 0..10..0. 

15.15.1 

1.. 0..3..0. 

16.17.1 

1.. 0.. 4.5.0. 

18.18.1 

3.5. 0. . 5. .0. 

19.20.1 

0.044,0., 0.056,0. 

21 . 21.1 

0.,0.,1.E15,0 

$ BLOCK R: 

END 
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36.. 36.. 8. .0.25.0.25. 20. 

386.3,5000. 

0.283,0.283,30000000. ,30000000. ,0.3, 0.3, 16. , 16 

0 . 001 , 20 . , 20 . , 8 . , 8 . 

620.. 1..1..1.. 1.25 

0.0001 3,0.000001 8,7 100., 120. 


Input Data Examples: Dynamic Load Minimization 

The input data for a design example of the dynamic load minimization problem is 
presented in Table 2. 


Table 2 


$ BLOCK A* 

OPTIMIZATION DESIGN OF DYNAMIC LOADING OF SPUR GEAR TOOTH 
$ BLOCK B: 

2,5 

$ BLOCK C: 

0,5,7,1000 
$ BLOCK D: 

$ BLOCK E: 

$ BLOCK F: 

5, 2, -1.0 

$ BLOCK G: 

30.. 50.. 36. 

8 .. 12 .. 8 . 

0.2,0.4,0.25 

$ BLOCK H: 

1.3 

1.4 

2.5 

3.6 

3.7 

$ BLOCK I: 

10 

$ BLOCK J: 

9.10.1 

0., 0., 180000. ,0. 

11 . 11.1 

0.,0.,1.E15,0. 

12 , 12,1 

-l.E+10,0.,0.,0. 

13,13,1 
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0.,0.,275.,0. 

14,14,1 

l.,0.,10.,0. 


15.15.1 

1.. 0..3..0. 

16.17.1 

1.. 0.. 4.5,0. 

18.18.1 

3.5.0. .5..0. 

19.20.1 

0.044,0., 0.056,0. 

21 . 21.1 


0.,0.,1.E15,0. 
$ BLOCK R: 
END 


36.. 36.. 8.. 0.25. 20. 

386.3,5000. 

0.283, 0.283, 30000000., 30000000., 0.3, 0.3, 16., 16. 

0 . 001 , 20 ., 20 ., 8 . , 8 . 

620.. 1..1..1.. 1.25 
0.00013, 0.0000018, 7100., 120. 


Program Operation 

After compiling, FORTRAN file SPUR is linked with COPES and ADS to obtain the 
completed spur gear optimization design program, where, file SPUR includes subroutine 
ANALIZ, TREAT, BENDST, GEOFAC, GEOFUN, GEODF, DYLOAD, SURDUR, 
SCORIN, FLASH, FILMSC, WEIGH, RATIOC, OUTPUT, LIFESG and FILMTH, i.e., 
LINK SPUR, COPES, ADS 

After linking, the input and output data file should be defined, then, the spur gear 
optimization design program could be run to obtain the optimal results. 

4.3 Design Examples 

Two optimization problems are considered here: life maximization and dynamic load 
minimization. 


Life Maximization 

The operation conditions are: 
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Speed (pinion RPM): 5000 
Normal load (LbS): 386.3 

Ambient temperature (°F): 120.0 


Gear material constants 

pinion 


gear 

Young’s modulus (psi) 

30 x 10 

6 

30 x : 

Poisson’s ratio 

0.3 


0.3 

Density (Lbm/in ) 

0.283 


0.283 

Lubrication constants 




Viscosity at 170°F (lb. sec/in 2 ) 


1.8 

x io' 6 

2 

Pressure-viscosity coefficient (in /lb) 

1.3 

o 

X 

Temperature-viscosity coefficient (deg R) 

7.1 

X 10 3 

Gear geometry constants 

pinion 


gear 

RMS surface roughness (micin) 

16.0 


16.0 

Pressure angle (deg) 

20.0 


20.0 

Measured error (in.) 

0.001 


0.001 

Length of shaft (in.) 

20.0 


20.0 

Diameter of shaft (in) 

8.0 


8.0 


Brinell hardness 
Life factor C L 
Hardness-ratio factor C. 
Reliability factor C R 
Temperature factor C T 


620.0 

1.0 

1.0 

1.25 

1.0 
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The life maximization problem is 


L(x) (Gear mesh life) 

1) 0<G (bending stress) <180000 

2) o H (Hertzian stress) < S H (fatigue strength) 

3) Scoring number < critical scoring number 

4) 0 < T f (Flash temperature) <275 

5) 1 < X (specific film thickness) <10 

6) 1.4 < nip (contact ratio) < 2.0 

7) 1.0 < W p (pinion weight) < 4.5 

8) 1.0 <W g (gear weight) <4.5 

9) 3.5 < C (center distance) < 5.0 

10) 0.044 < (length to diameter ratio) < 0.056 

11) Np t (for no involute interference) > 1 

12) Nj<N 2 

Side Constraints 

1) 30 <Xj, X 2 (the number of teeth) <60 

2) 8 < X 3 (diametral pitch) < 12. 

3) 0.2 < X 4 , X 5 (facewidth) < 0.4 

the initial design was 

Nj = N 2 = 36, P = 8, f 2 = f 2 = 0.25 (in.) 
has a mesh life of 666.1 million pinion rotations. 

The optimization algorithm used was the Modified Feasible directions method with 
polynomial interpolation with bounds for the one-dimensional search. After optimization, 
the final design has a mesh life of 1963.6 million pinion rotations. 

The final value of the design variables are 


Maximize 
Subject to 
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Nj = N 2 = 58.8 


P= 11.774 


f 1 = f 2 = 0.28 (in.) 


The initial and final dynamic loads, Wd, were 1141.1 and 1190.1, respectively. 

Dynamic Load Minimization 

For the dynamic load minimization problem, all the input constants, constraint 
conditions, side constraints and initial values of design variables are the same as in life 
maximization problem. Only the objective function becomes dynamic load Wd now. 

That is 

Minimize Wd 

the initial design has a dynamic load of 1 141.1 lb. 

After optimization, where the same optimization method as in life maximization 
problem is used, the final design has a dynamic load of 1089.2 lb. 

The final value of the design variables are 

Nj N 2 = 39.8 P = 8, fj = f 2 = 0.219 (in.) 


Summary 

It is seen from the results of these examples that the mesh life of spur gears can be 
increased significantly (195%) by increasing the number of teeth, diametral pitch and face 
width without any significant change in the weight and size of the gears. Meanwhile, all 
the constraints conditions considered here are satisfied. 

On the other hand, the dynamic load on spur gear tooth is reduced only a little (4.6%) 
by increasing the number of teeth and decreasing face width. The reason may be that the 
dynamic load on gear tooth mainly depends on the dynamic parameters such as speed and 
measured error on pair of mating teeth rather than the geometric parameters which are 
design variables here. 
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The nature of the dynamic load optimization problem needs to be investigated further. 
More reasonable design variables are expected to be chosen such that the dynamic load can 
be minimized much more. Besides, the constraint conditions corresponding to new design 
variables need to be developed. 

5.0 SPIRAL BEVEL GEAR ANALYSIS AND OPTIMIZATION 

Spiral bevel gears have advantages than other bevel gears. The mesh of spiral bevel 
gears has a rolling contact. Since spiral bevel gears have curved oblique teeth, the contact 
is smooth and continuous from end to end. Also as a result of their overlapping teeth 
action, spiral bevel gears will transmit motion more smoothly than spur gears and other 
bevel gears thereby reducing noise and vibration which become especially noticeable at 
high loads and high speeds. 

On the other hand, because of the complicated geometry of spiral bevel gear, the basic 
mechanisms which govern the major modes of spiral bevel gears are still not fully i 
understood. Thus the design and analysis of a spiral bevel gear set remains a difficult 
problem. 

The analysis presented is for the design of a set standard spiral bevel gear mesh. A 
typical objective is to maximize gear life without increasing gear weight. The transmitted 
horsepower, rpm of the pinion, and gear reduction ratio are given as input. A typical set of 
constraints includes bending fatigue, surface pitting, and scoring failure. Kinematic limits 
of face contact ratio, involute interference and other geometric limits are also considered. 

The study here is based on a thorough study of the kinematics of the gear mesh, such 
as those by Buckingham [9]. A Tregolds approximation to calculate spiral bevel gear 
geometry is used. The study is also based on a standard spiral bevel gear set with 
dimensions of the AGMA. 

The purpose is to establish an optimal design procedure for a set of standard spiral 
bevel gears which are given specified pinion speeds, input horsepower and gear reduction 
ratio. 
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Analytical Model for the Life of Spiral Bevel Gears 
Dynamic Load Bending Strength 

The basic geometrical parameters describing the spiral bevel gear are shown in figure 
5. The fundamental equation for bending stress in a gear tooth was developed by Wilfred 
Lewis in 1893. The AGMA has modified the basic bending stress by several factors that 
deal with the characteristics of specific application as follows: 

W t K„ P, KK 

s r~T'“ (86) 

K, 

s , sa » rr < 87 > 


where 

W ( = tangential load at the pitch diameter of large end 



126000 


% d P 


H p x 33000 
V 


(Lb) 


( 88 ) 


n p = Rpm of pinion (rpm 

d p = pitch diameter of pinion at large end (in) 

V = pitch line velocity (in/sec) 

Hp = input horsepower 

P d = diametral pitch (1/in) 

F = face width (in) 

J = geometric factor for bending 

K 0 = overload factor 

K y = dynamic factor 

K = size factor 
s 
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K = load distribution factor 
m 

K l = life factor 
K t = temperature factor 
K r = safety factor 

= calculated tensile stress at root of the teeth (psi) 
sat = allowable stress for material (psi) 


Usually sat = 30000 psi for an aircraft power drive case [10]. In general, gears should be 
designed for equal stress on gear and pinion. 

All of the K factors can be chosen in different cases by the user [10]. 

The geometry factor J is an index of the tooth geometry and stress concentration in 
the root fillet. It strongly depends on the cutting tool geometry. It is also an important 
factor in calculating the tensile stress. 

The AGMA standard has given the geometry factor J for a spiral bevel gear tooth at 
shaft angle 90‘, pressure angle 20°, spiral angle 35* graphically. In order to compute J for 
other cases, the polynomial approximation is used to obtain the J factor from the AGMA 
charts. 

A linear polynomial equation is used as follows: 


J - a 0 + a 1 N p + a 2^G + a 3 N P^G 


(89) 


In order to determine the coefficients a^a^a^a^ . . ., 5-10 points are chosen from each 
curve. Then, the least squares method is used to determine the a^a^ This gives 


J = 0. 09003 + 3. 9145 x 10 3 N n + 2. 544 x 10 3 N^ 

P U 

- 3. 2423 x 10 5 N p N G ( 9 °) 
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(12<N p <100 


12 < N G < 100) 


where N p = number of teeth on pinion 
Nq = number of teeth on gear 

This equation is not only convenient for computer code, but also is quite accurate for 
calculating geometry factor J. 


Pitting 

Another mode of failure is called "pitting" of the tooth surfaces. It is caused by too 
high contact stress known as "Hertz stress" at the contact surface. Here we assume that the 
magnitude of Hertz stress is a reasonable measure of the tendency of the surface to pit. 

In the AGMA standard, the Hertz stress can be calculated by the following formula 


s c = c p 


I 


w c c c c 

" t^o s m^f 

C„ d- F- I 


(91) 


C C 

S < s 

d c ~ a ac c C 
^TR 


(92) 


where 

W = tangential load at the pitch diameter of large end (Lb) 

C = elastic coefficient 
P 


C 


P 


1.5 




) 


p p = Poisson's ratio for pinion 
P G = Poisson's ratio for gear 
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Ep = Young's module for pinion 

E q = Young's module for gear 

d = pitch diameter for pinion at large end (in) 

F = face width (in) 

C Q = factor overload 

C = size factor 
s 

C = load distribution factor 
m 

C f = surface condition factor 
C y = dynamic factor 
= life factor 

C H = hardness ratio factor 
C T = temperature factor 
C k safety factor 
S p = calculated Hertz stress (psi) 

S = allowable Hertz stress (psi) 

= (0.4H b - 10) x 100 (psi) 

Hg = Brinell hardness of pinion 

H fi can be the natural hardness of the material or the limit of the tooth hardening. 

The geometry factor, I, depends on the radius of curvature at the point of contact and 
the load sharing between teeth. The AGMA standard has given the geometry factor I for a 
spiral bevel gear tooth at shaft angle 90°, pressure angle 20°, spiral angle 35* graphically. 
For the computer code, this number is obtained from the polynomial which fits the curve of 
the AGMA charts. The procedure is the same as for the geometry factor, J. The resulting 
polynomial equation is 
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I = 0. 0680 - 1. 0736 x 10 3 N p + 1 399 x 10 3 N Q 

- 1 235 x 10 6 N n N_ (93) 

P U 

(12 £ N p < 50 15 < N q < 100) 

All of the c factors can be chosen by user [5]. 

Scoring 

Besides breakage fatigue and pitting fatigue, the third mode of failure is scoring. 
Recently more attention has been paid to this phenomena. It usually occurs in high speed 
and high load conditions, with low viscosity lubricants when the overheating causes the oil 
film breakdown and the gear is kept in operation after pitting starts. This phenomena is 
also called metal-to-metal contact resulting in welding and tearing apart until the tooth 
surface is deteriorated. Actually, scoring is a lubrication failure, and is hard to analyze and 
predict. Research has reported that it is influenced by a combination of factors, such as 
pressure, sliding velocity, viscosity of the oil, temperature of the oil bath, properties of 
material and surface finish and treatment, etc. 

There are many ways to predict scoring risk. Two design approaches will be given 
here. The one efficient method is based on flash temperature (called hot scoring). The 
other is based on oil-film thickness (called cold scoring). 

The concept of flash temperature was first presented by Blok [16]as follows: 

T f = T fe + AT (94) 

Where: T fe is the temperature of gear blank in the contact zone. It is difficult to estimate. 
Often the blank temperature is approximated as the average of the oil temperature entering 
and leaving the gearbox [13]. In this paper, we choose T fe = 200°F. The increment of the 
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temperature, AT can be calculated as follows [13]: 


W 


t N 0. 75 


N 


AT = z t (-=rr ■ -p 


( 95 ) 


where 

W = tangential load of the pitch diameter at large end (Lb) 

F = face width (in) 

n p = rpm of pinion (Rpm) 

P d = diametral pitch (1/in) 

Z t = scoring geometry 


„ 0.0175(^-^TMG)P 4 d 

3/4 1/4 

(cos<i >t r 4 [p 1 p 2 /(p 1 +p 2 )] 1/4 


(use absolute value of Z { ) 
where 

<t> t = transverse presure angle tan<J> t = tan<t)/cosv 

\)/ = spiral angle 

<J> = normal pressure angle 

MG = gearratio 

p l = the curvature of contact point for pinion 

p 2 the curvature of contact point for gear 

(Note: Here we use the material that is steel on steel and take the coefficient of friction as 
0 . 06 ). 

Because scoring is considered most likely to occur at the lowest point of single tooth 
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contact on the pinion or the highest point of single tooth contact on the pinion, these two 
points are checked. Normally, scoring calculations are made for the pinion only and the 
gear is handled indirectly. If the pinion is safe the gear should be safe. 

The formula used here to calculate the lower single tooth contact curvature for spur 
and helical gears at lowest single tooth contact is, 


p l = J T l~ ( rcos< t> t ) 2 - pcos^t (97) 


p 2 = cos<|> t - Pj 


(98) 


at highest single tooth contact: 


pj cos<(> t = p 2 


(99) 


p 2 



(R cos <j) t ) 2 


p cos <|) 


( 100 ) 


Where 

r Q = outer radius for pinion (in) 

r = pitch radius for pinion (in) 

R q = outer radius for gear (in) 

r = pitch radius for gear (in) 

p = circular pitch p = 7t/P d 

C = center distance (in) 

P = transverse pressure angle 

Because the spiral bevel gear has a complicated tooth geometry which is difficult to 
deal with analytically, in order to calculate the two critical point curvatures, a Tregold’s 
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approximation which reduces the spiral bevel gear problem to the helical gear is used here 
[9]. Thus we can use the helical gear formula to solve the spiral bevel gear problem. It 
uses the equivalent helical gears of the back cones of the bevel gears for tooth numbers 
greater than 8. This method is sufficiently accurate for practical purposes. The equivalent 
helical gear is calculated as [9]: 

R vp = cos y A R vg = cos y B (101) 

= pitch radius of equivalent helical pinion (in) 

R yc = pitch radius of equivalent helical gear (in) 

r = pitch radius of bevel pinion at the large end (in) 

R = pitch radius of bevel gear at the large end (in) 

y = pitch angle for pinion 
A 

y = pitch angle for gear 
fi 

Since bevel gears do not have any center distance dimension, the equations (l)-(4) can be 
modified as follows: 

At lowest single contact point 

Pi = V^o-CRvp^t) 2 “ Pcos<|> t (102) 


P 2 = (R vp + R VG )cos ^ t 


(103) 


At highest single contact point 
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Pl = ( R vp +R V g> C 0 S< t > t - P2 


(104) 


P 2 



p cos <]> 


(105) 


where 


r 


0 



1-0/Pd 


(106) 


R=R V „ + 1.0/P. 


k 0 


vg 


(107) 


Tq = outer radius of equivalent helical pinion 
R q = outer radius of equivalent helical gear 

Now we can calculate Z t at the two critical points. Then choose the smaller Z { to calculate 
AT to get the critical flash temperature. 

For aircraft spiral bevel gears, an allowable surface flash temperature is 493*F [17]. 
Here 400°F is imposed as an upper bound limit. 

The other method to judge the risk of scoring is the scoring criterion number. The 
formula is: 


W. 


" t 3/4 n P 
Calculate scoring number = Cp— ) — j- 


The AGMA standard has published the scoring criterion number table for various oils at 
various gear blank temperatures. If the calculated scoring number is above the value 
shown in the table, there is a risk of scoring [13]. 

Cold scoring occurs when the elastohydrodynjamic oil film thickness is small 
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compared with the surface roughness and the lubricant does not have enough additives to 
prevent scoring on the contact tooth surfaces. A criterion used to determine the possibility 
of surface distress is the ratio of film thickness to composite surface roughness [13]: 


X = 


min 

a 


>1 


(108) 


where ° = V°P + *« 

X = film parameter 

h . = minimum oil film thickness (in) 

<7 p = surface roughness of pinion (rms) 

a = surface roughness of gear (rms) 

s 

There is a risk of cold scoring if X is less than 1. Here the lower bound of 1 is used. 

The calculation of h min is a quite complicated problem. It requires special data about 
the oil. A simple but approximate calculation for the minimum oil film thickness at the 
pitch line is used here 

44. 6 r e (lubricant factor) (velocity factor) 

^min ” (loading factor) (M^ n ) (109) 


where 

r e = effective radius of curvature at pitch diameter 


r 


e 


csin( l>n MG 
cos 2 \]/ (MG + 1) 2 


(HO) 
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(Ill) 


Spiral Bevel gears do not have any contact distance, so use c = (MG + l)r 

r sin K MG 

e COS 2 y ^G 1 


r = pitch diameter at large end for pinion 

<)> n = normal pressure angle 

\| / = spiral angle 

MG = gearratio 

Lubricant factor = (aE') 0,54 

a = lubricant pressure-viscosity coefficient in /lb 

E ' = effective elastic modulus for a steel gear set 

= 4 = 51. 7 x 10 6 psi . 

2 ( 1 — v 2 ) 


E = 3 x 10 psi 


v= 0.3 


Ho‘ 


0.70 


Velocity factor = ( g. ) 


(113) 


|i o = lubricant viscosity at operating temperature (cp) and u = rolling velocity at the pitch 
line 


K n p d sin (|) t 

u = in / second 


^0 H*e[|3+ ( t+46 o t +460 )] 


(114) 
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|i = viscosity at temperature T Q 

P = lube, temp- viscosity coefficient at lube ambient temperature (deg. R) 

T 0 = tooth bulk temperature 

T = temperature (deg. F) 

d = pinion pitch diameter at large end (in) 

<J> t = transverse pressure angle 

Loading factor = (Wj/FE'r^ 0,13 
W = tangential load (Lb) 

F = face width (in) 

Face Contact Ratio 

In order to insure a quiet, smooth, reliable gear mesh, the face contact ratio must be 
considered. The face contact ratio for a given spiral bevel gear is the ratio of face advance 
to the circular pitch. It must be greater than a minimum number. The value 1.25 is 
commonly used for spiral bevel gear. 

Involute Intereference 

Involute interference occurs when the addendum circle of one gear crosses the line of 
action past its point of tangency with the base circle on its mating gear. That is, its contact 
occurs below the base circle of the pinion on the non-involute portion. This phenomena 
also is called "undercut" and this weakens the gear. There are many ways to avoid 
undercut. The use of minimum pinion teeth to eliminate interference is imposed here. 

The lower bound of teeth number for pinion is 12 for spiral bevel gears, when the 
pressure angle is 20*. On the other side, in order to avoid undercut with low numbers of 
teeth and balance the strength of gear and pinion teeth, there is a limit for gear and pinion 
teeth: 
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N + N_ > 34 
P G 


(for pressure angle equals 20 °) 


(115) 


Ratio of Face Width to Cone Distance 

The AGMA standard recommended the ratio of face width to cone distance is equal to 
V 3 or P yiO, whichever is smaller. 

5.1 The Spiral Bevel Gear Program 

The spiral bevel gear analysis program is written in the same general format as the 
spur gear analysis program. That is, the main routine is a sub routine called ANALIZ 
which can be directly coupled with the COPES/ADS optimization program. 

The analysis routines include ANALIZ, INPUT, OUTPUT, PITTING, BENDING, 
SCORING, CONTR, WEIGHT, FLSHT, FMTH, LFAFBG, and LIFEBG. These are all 
contained in a FORTRAN file called SPIRLB. Subroutine ANALIZ is the main analysis 
control routine which is called by COPES for optimization. This routine contains the 
GLOBCM common block which transfers data between the analysis and optimization. 
Subroutine INPUT reads all analysis data and subroutine OUTPUT prints the analysis 
results. 

5.2 Parameters Contained in the Global Common Block 

In subroutine ANALIZ, the labeled common block, GLOBCM, contains the 
information that is transferred between the analysis and optimization routines. The 
common block is as follows: 

COMMON/GLOBCM/NTA,NTB,FWAD,FWBD, PR, SIGMA, PSI, PHI, SC,ST, SCORN, 
FLASHT,LAMDA,FCON,MF, WEIGHT, INTF1,INTF2,LF, MG 

These parameters are defined in the following table for convenient reference: 
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Fortran 


Loc Parameter Name Units Definition 


1 

n p 

NTA 


Number of teeth in pinion 

2 

N o 

NTB 


Number of teeth in gear 

3 

F 

FWAD 

in. 

Face width of pinion 

4 

F 

FWBD 

in. 

Face width of gear 

5 

p d 

PR 

l/in. 

Diametral pitch 

6 


SIGMA 

degrees 

Shaft angle 

7 


PSI 

degrees 

Spiral angle 

S 


PHI 

degrees 

Pressure angle 

9 

S 

c 

SC 

psi 

Hertzian stress 

10 

s , 

ST 

psi 

Bending stress 

11 

SCON 

SCON 


Scoring number 

12 

T 

f 

FLASHT 

deg. F 

Rash temperature 

13 


LAMDA 


Thickness/roughness ratio 

14 

FCON 

FCON 


Ratio of face width/cone dist. 

15 

Mp 

MF 


Face contact ratio 

16 

WEIGHT 

WEIGHT 

Imb 

Gear set weight 

17 

INF1 

INF1 


No involute inter, lower bound 

18 

INF2 

INF2 


Sum of gear mesh 

19 

l f 

LF 

10 6 rev. 

Gear mesh life 

20 

m g 

MG 


Gear ratio 


Input Data Description 

As with the spur gear design, the data file consists of two parts. The first is the 
COPES input data which defines the optimization problem and the second is the analysis 
data for the spiral bevel gear. 
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Spiral Bevel Gear Input Data 

There are seven lines of data as defined below. 

Line 1 

GEMATA, GEMATB, LUBNA 
GEMATA - pinion material 
GEMATB - gear material name 


LUBNA 

- lubricant name 

Lind 

RPMA,SC,ST,Y1ESA,Y1ESB,YMAD,YMBD 

RPMA 

- RPM of pinion 

SC 

- allowable contact stress (psi) 

ST 

- allowable bending stress (psi) 

Y1ESA 

- pinion material yield stress (psi) 

Y1ESB 

- gear material yield stress (psi) 

YMAD 

- pinion material Young’s modulus (psi) 

YMBD 

- gear material Young’s modulus (psi) 

Ung 2 

NTA,MG,HP,HB 

,SRMS1,SRMS2,BETA,TEMLUB,SPHEA,SPHEB 

NTA 

- number of teeth in the pinion 

MG 

- gearratio 

HP 

- input horsepower 

HB 

- minimumb rinell hardness of the gear set 

SRMS1 

- pinion rms surface roughness (micro-inches) 

SRMS2 

- gear rms surface roughness (micro-inches) 

BETA 

- lube temp-viscosity coef. at lub ambient temp (deg. R) 

TEMLUB 

- lube ambient temperature (deg. F) 
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SPHEA - specific heat of the pinion 

SPHEB - specific heat of the gear 


Line 4 

ALPHA, VlSD,DENSA,DENSB,CONDTA,CONDTB 


ALPHA 

2 

- lub pressure-velocity coefficient (in /lb) 

V1SD 

- lube viscosity at lube ambient temperature 

DENSA 

- pinion mass density (lbm/in ) 

DENSB 

3 

- gear mass density (lbm/in ) 

CONDTA 

- pinion thermal conductivity 

CONDTB 

- gear thermal conductivity 

Line 5 

FW AD, FWBD,PR, SIGMA, PSI, PHI 

FWAD 

- pinion face width 

FWBD 

- gear face width 

PR 

- diametral pitch 

SIGMA 

- shaft angle 

PSI 

- spiral angle 

PHI 

- pressure angle 

Line 6 

CA,CS ,CM,DF,CL,CH,CT,CR 

CA 

- overload factor 

CS 

- size factor 

CM 

- load distribution factor 

CF 

- surface condition factor 

CL 

- life factor 
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CH 

- hardness ratio factor 

cr 

- temperature factor 

CR 

- safety factor 

Line 7 

KO,KM,KL,KT, 

KR,FRM,EPSA,EPSB 

KO 

- overall factor 

KM 

- load distribution factor 

KL 

- life factor 

KT 

- temperature factor 

KR 

- safety factor 

FRM 

- friction coefficient of material 

EPSA 

- Poisson’s ratio for the pinion 

EPSB 

- Poisson's ratio for the gear 


5.3 Design Examples 

The typical spiral bevel gear optimization problem can now be stated as 


Maximize 

Life 

Subject to 

(1) 

Constraints to 

(2) 


(3) 


(4) 


(5) 


(6) 


(7) 1 


(8) ' 


(9) 


(10) Sum of number for pinion and gear 
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In the following example the initial design was assumed that a single reduction spiral bevel 
gear set is to be designed for a 450 horsepower input from pinion shaft, the RPM of the 
pinion is to be 15000 RPM. 

The material AISI 9310 is chosen here. The gears are assumed to be made from case 
hardened alloy steel with ground tooth surfaces. Its hardness is 620 H fi . A constant 
pinion speed and gear ratio are maintained during the optimization. 

Super-refined naphthenic mineral oil is used as the lubricant. 

The initial values are described as follows: 


Gear Data: 


Number of teeth for pinion 

31 

Number of teeth for gear 

46.5 

Gearratio 

1.5 

Shaft angle 

90 degrees 

Pressure angle 

20 degrees 

Spiral angle 

35 degrees 

Diametral pitch 

6/in 

Standard Operating Conditions 


Pinion 

1500 Rpm 

Gear 

1000 Rpm 

Ambient temperature 

100°F 

Gear Material Data 


Material 

AISI 93 10 

Density 

0.283 Lb/in 3 

Young's modulus 

30,000,000 psi 

Poisson's ratio 

0.3 
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.6 RMS 


Surface finish 
Lubricant Data : 

Materials 

Temperature-viscosity coefficient 
Thermal conductivity at 100°F 
Specific heat 

Pressure viscosity coefficient 
Lub viscosity at ambient temperature 


Super-refined, naphthenic, mineral oil 
7100°K 

5.84 cp (centipoint) 

956 

0.00013 in 1 2 /Lb 
64.7 cp 


Here two examples are presented. Both examples have the same design variables, 
constraints and objective function. The only difference is in the bounds on the gear weight, 
and the initial values of the design variables. Also, example 2 is for a different gear ratio. 

The number of pinion teeth, face width for (same face width for gear) and diametral 
pitch are design variables for both examples. 

The objective function is to maximize the gear mesh life. 

Ten sets of constraints, which define the design space, are used in both examples. 
Their lower and upper bounds are as follows: 

1 . 0 <, Bending stress < 240000 Lb/in 2 . 

2. 0 < Pitting stress < 180000 Lb/in 2 5 . 

3. 0 £ Scoring number < 17000. 

4. 0 < Flash temperature < 400. 

5. 1 < Film parameter < 10 6 7 8 9 10 . 

6. 0 < Ratio of face width to cone distance <10. 

7. 1.25 < Face contact ratio <10. 

8. 16 < Total weight < 17.30 for example 1. 

13 < Total weight < 13.27 for example 2. 

9. 12 < Involute interference < 50. 

10. 40 < Sum of number pinion and gear teeth < 150. 
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The optimization results for example 1 are listed in Table 1, where the initial number 
of pinion teeth is 31. The gear ratio is 1.5 so the number of gear teeth is 46.5. The initial 
value for face width is 2.5 in, the initial diametral pitch is 6/in. In Tables 1 and 2, the 
parameters ISTRAT, IOPT and IONED refer to the strategy, optimizer and one- 
dimensional search options in the ADS program [5]. 

The initial design is an infeasible design. The upper bound of ratio of face width to 
cone distance is violated. The initial value for gear mesh life is 12958 x 10 6 cycles. After 
optimization, the optimizer brought the initial infeasible design to a feasible design. Here 
four different combinations of strategies have been used on this example and none of the 
constraints were violated. The gear mesh life is increased. Strategy 8 (Sequential 
Quadratic Programming) failed to produce an optimum design. 

In example 2, different initial values of design variables were used. The number of 
the pinion teeth is 18, and the gear ratio is 3.5. The optimization results for example 2 are 
shown in Table 2. It is seen that the initial design is also infeasible. Two constraints are 
violated, one is the upper bound of flash temperature, the other is the upper bound on the 
ratio of the face width to cone distance. The optimizer brought the initial infeasible design 
to a feasible optimum design. At the optimum, none of the constraints were violated. The 
life was reduced from its initial infeasible value, but is optimum for the constraints that 
were imposed. Strategy 8 again failed to produce the optimum, although it came close. 
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Table 1 
EXAMPLE 1 


ISTRAT 


0 

6 

8 

6 

IOPT 


5 

5 

5 

4 

IONED 


7 

7 

7 

5 

IPRINT 


3050 

3500 

3500 

3500 


INITIAL 


OPTIMUM 



n p 

31 

16 

16 

22.37 

16 

n g 

46.5 

24 

24 

33.55 

24 

FW 

2.5 

3.131 

3.15 

2.331 

3.146 

PR 

6 

3.17 

3.18 

4.29 

3.18 

MG 

1.5 

1.5 

1.5 

1.5 

1.5 

Weight 

17.25 

17.30 

17.22 

17.225 

17.25 

Life 

12958 

20695 

20562 

13547 

20501 

Active 

0 

3 

1 

0 

3 

Violated 

1 

0 

0 

0 

0 
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Table 2 
EXAMPLE 2 


ISTRAT 


0 

6 

8 

6 

IOPT 


5 

5 

5 

4 

IONED 


7 

7 

7 

5 

IPRINT 

INITIAL 

3050 

3500 

OPTIMUM 

3500 

3500 

n p 

18 

12 

12 

12 

12 

n g 


42 

42 

42 

42 

FW 

2.5 

2.504 

2.499 

2.501 

2.501 

PR 

6 

4 

4 

4 

4 

MG 

3.5 

3.5 

3.5 

3.5 

3.5 

Weight 

13.24 

13.26 

13.209 

13.239 

13.259 

Life 

29658 

26908 

26733 25668 

26847 

Active 

0 

3 

3 

3 

3 

Violated 

2 

0 

0 

0 

0 
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Summary 

From the results of the examples that are presented here, it can be seen that the life 
can be increased by decreasing the diametral pitch without increasing the weight. Also the 
face width is another significant factor for life. The optimizer changes the design variables 
until an acceptable design can be gained. In these two examples, the optimizer brought the 
initial infeasible design to feasible (even though sometimes the life is reduced). In example 
1, after optimization the final mesh life is roughly 59% more than the initial design. In 
example 2 after optimization the final mesh life is roughly 9% less than the initial design, 
but it is feasible design. 

Further Work 

This program is only used for a standard bevel gear set with shaft angle 90*, pressure 
angle 20* and spiral angle 35* because of lack of some information. It may be expected that 
the spiral angle also is a significant factor for the mesh life. Thus, if enough information is 
available to modify the program and include this effect, the use of this program will be 
broadened. 


6.0 SUMMARY 

Design of precision gears has been addressed as a numerical optimization task. The 
principal objective here has been to investigate how optimization methods may be applied to 
this design problem and to demonstrate the resulting technology. To achieve this, the 
development of software that demonstrates the theoretical concepts is necessary. This is 
consistent with the overall goal of research, which is to produce a usable tool for 
engineering design. 

Two gear design programs have been developed, one for spur gear design and one 
for spiral bevel gear design. These codes have been demonstrated with test cases to 
demonstrate their applicability to typical gear design problems. The key to successful 
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design optimization is that a wide variety of constraints must be considered. For example, 
if a gear is designed for minimum noise without regard to bending or Hertzian stress 
constraints, the resulting design will be of little practical use, since it will soon fail due to 
insufficient strength. Thus, an essential part of this effort has been to consider many 
constraints, something that is natural for the numerical optimization approach to design. 

A basic gear design capability has been developed here. This effort has dealt only 
with the gear itself, as opposed to the overall speed reduction system which includes 
shafting, housing, powerplant, etc. Also, the analysis tools used here are relatively basic. 
Continuing research is considering more detail in the analysis. This includes more 
sophisticated models for noise generation as well as more detailed models of the gear teeth 
themselves. In the later case, research is being undertaking to model the gear teeth using 
finite element methods. An example of this effort is to design gear teeth which have an 
involute profile under the deformed state , as opposed to the design where the unloaded gear 
is an involute shape. 

The use of numerical optimization has been shown to significantly improve the life of 
gears, reduce noise, and otherwise enhance the design process. Most importantly, this 
technology provides the ability to investigate the tradeoffs between the competing 
constraints in a rational way. 
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FIG. 1 - Spur Gear Geometry 



FIG. 2 - Geometric Details 
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FIG. 5 - Spiral Bevel Gear Geometry 
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